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Abstract. For first order systems, we obtain an efficient bound on the ex- 
ponential decay of an eigenfunction in terms of the distance between the cor- 
responding eigenvalue and the essential spectrum. As an example, the Dirac 
operator is considered. 



1. Introduction 

There are two different methods for proving exponential decay of eigenfunctions 
ip for second order selfadjoint elliptic operators H . The first of them is due to Shnol' 
[Sj who proved that an eigenfunction corresponding to an eigenvalue A satisfies the 
estimate 

/ \iP{x)\ 2 e 25lxl dx < oo (1.1) 

where 8 depends on the distance d(X) = dist{A, a eS s{H)} between this eigenvalue 
and the essential spectrum a ess (H) of the operator H . More precisely, it was shown 
in [5] that, for the Schrodinger operator H — — A + V(x) with a semibounded real 
potential in the space L2(M. d ), estimate (|1.1[> holds with any S such that 

S < 2- 1 ln(l + d 2 (A)a- 2 ) 

where a depends only on 

Vq = inf V(x). 

xen d 

Another method originated in papers [H [5] where eigenfunctions corresponding 
to eigenvalues lying below <7 ess (H) were considered for the iV-particle Schrodinger 
operator H . In this case estimate (11.11) with any 5 < -J d(X) was proven. Then 
Agmon [I] extended this result to very general differential operators of second order. 
This estimate is of course more efficient and stronger than that of Shnol' but, in 
general, does not hold for eigenvalues lying in gaps of a ess (H). This can be shown 
on an example of an one-dimensional Schrodinger operator whose potential is a sum 
of a periodic function and of a function with compact support. From a technical 
point of view, an estimate of eigenfunctions for eigenvalues in gaps is essentially 
more difficult than for eigenvalues below a ess (H) because in the first case one has 
to work with the expression \\Hf — A/|| whereas it suffices to consider the quadratic 
form (Hf, /) in the second case. 

We note also paper [3] where an estimate of eigenfunctions for eigenvalues in 
gaps was obtained for random Schrodinger operators. 
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Our goal here is to study exponential decay of eigenfunctions for first order 
matrix differential operators 

H =- i Y, A ^ + V ( x ) ( L2 ) 

j=l 3 

acting in the space H = L2(M d ;C"). Here Aj — A*, j — 1, . . . , d, are constant 
matrices and V(x) is a symmetric matrix-valued functiion. Set 

d 

J = TD3x\y2A j X j \, || = He- ( L3 ) 

1x1=1 u 

For example, 7=1 for the Dirac operator. Our main result is the estimate 

/ |^(s)| 2 e M l s ldx < 00, \/S<j- 1 d{\), (1.4) 

for all eigenvalues (including those lying in gaps of a ess {H)). 

Of course for first order differential operators the problem considered is essen- 
tially simpler than for operators of second order. This allows us to follow sufficiently 
closely the method of pQ and to obtain estimate (|1.4j) without any specific assump- 
tions on the coefficients Aj and V(x). Note that for second order elliptic operators 
one can pass from weighted /^-estimate (jl.ll) to the pointwise estimate 

\ip(x,\)\<Ce- 5 W, 8<^/d(X}. 

Such a passage is eventually possible also for first order systems but would require 
more or less stringent assumptions on coefficients of the system. 

Here and below C denotes constants whose precise values are irrelevant; / is the 
identity operator (in different spaces). 



2. The main result 

Let us first discuss the precise definition of the operator H . Derivatives in (|1.2p 
are understood in the sense of distributions. We suppose that V(x) is a locally 
bounded function, but we do not make any specific assumptions about its behavior 
at infinity. Instead of this we a priori assume that H is selfadjoint on a domain 
with the following natural properties. 

Assumption 2.1. Operator (| 1 . 2[) is selfadjoint on domain T>(H) such that 
1° If a (vector-valued) function f(x) has compact support and 

d 

^Ajdf/dxj e i 2 (K d ;C"), 

3=1 

then / G V{H). 

2° If / G V(H), then 

d 

^Ajdf/dxj g 4' oc) (^ d ;C"). 

i=i 
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If the function V(x) is bounded, then the operator H is sclfadjoint on the set of 
functions / £ Ti such that 

d 



/ | Ajdf I dxj | 2 dx < oo. 



In this case Assumption 12.11 is of course satisfied. Recall that the operator H is 
called elliptic if the matrix 



d 



A(0 = > A& (2.1) 



does not have zero eigenvalues for £ ^ or, equivalently, 

\A(Op\>c\ P \\Z\, oo, peC", (eK d . 

Elliptic differential operators with bounded V(x) are selfadjoint on the Sobolev 
class H 1 (]R d ;C n ) =: V(H). We emphasize however that we do not assume neither 
ellipticity of the operator H nor boundedness of V(x). 

The following standard result is a direct consequence of the spectral theorem. 

Lemma 2.2. For every e > there exists R such that 

||ff/-A/||>(<f(A)-e)||/||, d{X)=dist{X,a ess (H)}, (2.2) 
for all f £ 'D(H) such that f(x) = in the ball \x\ < R. 

Proof. Suppose to the contrary. Then for some e > and all R there exists f R £ 
V(H), \\f R \\ = 1, such that f R (x) = for |at| < R and 

\\Hf R - Xf R \\ < d{\) - e. (2.3) 

Let us choose some a < d(\) and set X = [X — a, X + a], g R = f R — E{X)f R where 
E(X) is the spectral projection of the operator H corresponding to the interval X. 
Since g R = E(M. \ X)g R , the spectral theorem implies that 

\\Hg R -Xg R \\>a\\g R \\. (2.4) 

Using that f R — > weakly and that the operator E{X) has finite rank, we find that 

11^-/^11=0(1), \\Hg R - Hf R \\ =o(l) 

as R — > oo. Thus, comparing estimates p.3| and (|2.4p . we arrive at the contradic- 
tion if a > d(X) — e and R is sufficiently large. □ 

The following theorem is our main result. 



Theorem 2.3. Let the operator H be defined by formula (|1.2[) . let Assumvtion Yl.W 
be satsified, and let 7 be constant (| 1 - 3(1 - // Hip — Xtp where tp £ T>(H), then 
estimate (11.411 holds. 



Proof. Put e = 2 1 (d( X) — jS) and choose R in such a way that estimate (|2.2p holds 
if / € £>(#) and /(x) = for |x| < R. Fix some number i?i > i?. Let C £ C°°(R + ), 
C(t") = for r < R and C(r-) = 1 for r > R\. For an arbitrary (large) R2 > Ri, we 
choose p £ C°°(E + ) such that p{r) = 6r for r < R 2 , p'{r) < 6 for r £ (R 2l 2R 2 ) and 
p(r) = const for r > 2i?2. Set n(r) — Q{r)e p ^ and /(x) = n(\x\)ip(x). Of course, 
f = ip — (1— i])ip £ T>(H) because 1 — 77 € C^°. It follows from equation Hip = Xpj 
that 

(ff/)(ar) - Xf(x) = -ir- 1 V(r)A(a;)V(a ; ) 
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where the matrix A(x) is defined by formula (12. Using notation (II. 3| . we find 
that 

||#/ - A/|| < 7 ||r?VII < 7(l|p'e P C^II + l|e p CVII)- (2.5) 
Since p'(r) < *, we have that 

\\p'e p Cn<6\\f\\. 
Since supp£' c [R, Rt], we have that 

||p/CVII < c\\4>\\ 

where the constant C does not depend on R 2 . Thus, it follows from (|2.5p that 

||if/-A/||< 7*||/|| +CU\\. (2.6) 

On the other hand, the function / satisfies estimate (|2.2j) because f(x) — for 
\x\ < R. Comparing estimates (|2.2I) where £ = 2~ 1 (<i(A) — 7*) and (|2.6[) . we see 
that 

(d(A)- 7 *)||/||<C||V|| 

and hence 

f \ip(x)\ 2 e 2S Wdx < C\\ipf. 

JR 1 <\x\<R 2 

This proves (|1.4p because C does not depend on R 2 which can be chosen arbitrary 
large. □ 

3. The Dirac operator 



As an example to which Theorem 12 .31 directly applies, we now consider the Dirac 
operator describing a relativistic particle of spin 1/2. Let Ti = L2(M 3 ; C 4 ) and 

3 g 

H = -i^2aj- \-ma + W{x), (3.1) 

OXj 

3=1 J 

where m > is the mass of a particle and 4x4- Dirac matrices a.j = a* satisfy 
the anticommutation relations 

ctjOtk +a k a j ; = 0, j^k, a? = I, j,k = 0,1,2,3. 

These relations determine the matrices otj up to a unitary equivalence in the space 
C 4 . For example, we can express them in terms of the Pauli matrices 

vi = ( I I ) , °* = ( ° 7 ) ' ff3 = ( J -1 ) (3 - 2) 

by the relation 

«*=(° 0)' z = 1 ' 2 ' 3 ' a °=(o -j)' (3 ' 3) 

We do not require any special assumptions on the matrix- valued function W(x). 
Operator (|3.1|) has form (|1.2|) where d = 3, n = 4, Aj = atj and 

V(x) — ma + W(x). 

It follows from formulas (|3.2|) . (13. 3j) that matrix (|2.1j) now equals 

A(£) -( ° a ^ 3 \ 
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Since (X)j=i a iij) 2 — I£| 2 j an easy calculation shows that the matrix A(£) has 
two eigenvalues |£| and — 1£| (of multiplicity 2 each) and hence expression (jl.3|) 
equals 1. Thus, according to Theorem 12.31 estimate (|1.4[) where 7 = 1 is true 
for all eigenfunctions of the Dirac operator. Perhaps this result is new even for 
perturbations W(x) decaying at infinity when a ess (H) = (—00, —to 2 ] U [to 2 , 00) has 
the gap (—to 2 , to 2 ). 

Recall that estimate (II. lj) with S < yj d(X) is sharp for the Schrodinger operator. 
We cannot claim the same about our estimate for the Dirac operator. To clarify 
these statements, let us consider the simplest one-dimensional case. 

Let first H — —d 2 /dx 2 + V(x) where V(x) is a function with compact support. If 
A < is an eigenvalue of H , then, for large ±x > 0, the corresponding eigenfunction 

equals ip(x) = c±e - V / I^Tl :E l where c± are some constants. Since now d(X) = |A|, one 
cannot hope to prove estimate (jl.lj) even for S = y/d(Xj. 

The one-dimensional Dirac operator is defined by the formula 

*=-<(; ;)<*/*»+(" _° m )+w 

in the space TC = £2(8; C 2 ). If the (matrix-valued) function W(x) has compact 
support and A € (—to 2 , to 2 ) is an eigenvalue of the operator H, then, for large \x\, 
the corresponding eigenfunction equals 

ip(x) = c ± (V^T\jV : ^~Xsgnx) t e- Vm2 - x2lxl , ±x > 0. 

Thus, the integral in (|1.4[) is convergent for all S < \J to 2 — A 2 which is better than 
5 < d(X) = m — |A| as required in (|1.4I) . 

4. The Hill operator 

First we discuss a curious application of the Agmon estimate to the one-dimensional 
Schrodinger operator 

H = -d 2 /dx 2 + V(x) (4.1) 

with a periodic potential V(x) — V(x + 1). Recall that for any A in the gap of the 
spectrum a(H) — o ess (H) of H the equation 

- y"{x) + V(x)y(x) = Xy(x) (4.2) 

has two non-trivial solutions y±(x) admitting the representation 

y±(x) = p* x p±(x) (4.3) 

where both functions p±(x) are either periodic or antiperiodic (p±(x) = p±(x + 1) 
ot p±(x) — — p±(x+l)) and the number p = p{X) > 1 is known as the multiplicator. 
Thus, the solution y±{x) exponentially decays as x — > ±00. 

Proposition 4.1. Suppose that 

X < A := inf a ess (H). 

Then 

In 2 p{X) > A - A. (4.4) 
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Proof. Let us consider the auxiliary family of operators 

H a = ~d 2 /dx 2 + V[x) - aQ{x) (4.5) 

where Q 6 Cg°(R) and Q(x) > 0. If Q(x) is not identically zero, then we can 
choose a > such that A is an eigenvalue of the operator H a . The corresponding 
eigenfunction 

ip(x) = c±y±(x) (4.6) 

for sufficiently large ±x > 0. On the other hand, it satisfies estimate (jl.lj) where 
5 is any number smaller than -\/Ao — A. By virtue of (|4.3p this is possible only if 
estimate (14.41) is true. □ 



Finally, we show that estimate (JTTTJ) with arbitrary 8 < \/d(X) is not in general 
true for eigenvalues in gaps. To that end, we need some additional, although always 
quite elementary, information on the Hill operator. Let (p{x, A) and 9(x, A) be the 
solutions of equation (|4.2p satisfying the conditions 

V5(0,A) = 0, </?'(0,A) = l, 6»(0, A) = 1, 0'(O,A) = O, 

and let 

^(A) = 2- 1 (/(1,A)+0(1,A)). 
Then |-F(A)| > 1 for regular points of H, and the equation |F(A)| = 1 determines 
band edges of its spectrum. The functions p(X) and F(X) are related by the formula 



p(\) = \F(\)\ + y/F*(\)-l. (4.7) 

We note that the functions p±(x) in (|4.3[) are periodic (antiperiodic) if F(X) > 1 
(ifF(A)<-l). 

Let A be a regular point close to one of two (one for A < cr(H)) band edges. 
Denote this edge by A*. It follows from (|4.7p that, as A — > A*, 



p(A) = 1 + y/2\F'(\.)(\ - X,)\ + 0(A - A,) 

and therefore 

In 2 p(X) = 2|F'(A*)| |A - A*| + 0(\X - A,| 3 / 2 ). (4.8) 
Let us now use standard asymptotic formulas 

<p'(l,A) = cos \/A + 0( -=), 9(1, X) =cosVx + o(-f=), X^oo, 



which can easily be deduced from the Volterra integral equations for solutions 
ip(x,X) and 6(x, A) of differential equation (|4.2p . Moreover, these asymptotic for- 
mula can be differentiated so that 

F , (A) = _sin^ + /ly A _ >Qo< 



2y/X U 

Comparing formulas (|4.8p and (|4.9|) . we obtain the following result 



Proposition 4.2. For an arbitrary c > 1, sufficiently distant gaps and X sufficiently 
close to a band edge A* the estimate 

ln 2 p(A) < cA7 1/2 |A - A„| 

is true. 

Let us again consider operator (|4.5p . We use the following simple result. 
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Lemma 4.3. For an arbitrary regular point A € M of the operator H , there exists 
Q = G 2 € Co°(K) and a£l such that the operator H a has the eigenvalue A. 

Proof. It suffices to show that, for suitable G and a, the Birman-Schwinger equation 

u- aG(H - \I)- 1 Gu = 0, H = H , (4.10) 

has a non-trivial solution u. Let us check first that we can find G such that the 
operator G(H — \I)~ 1 G ^ 0. Supposing to the contrary, we denote by G n the 
operator of multiplication by the function G n (x) = g(x/n) where g G Cq°(1R) and 
g(x) = 1 for |x| < 1. By our assumption, G n (H — AJ) _1 G n = and hence 

(G n (H - XI)- 1 G n u l ,u 2 ) = 

for all n and all Ui,tia € L20R). Since G„ converges strongly to the identity 
operator, we have that ((H — XI)~ 1 u\,U2) = 0. Since 1*2 is arbitrary, this implies 
that (H — A/) _1 iti = and hence u\ = Q although u\ is arbitrary. 

Let us pick an arbitrary eigenvalue \i 7^ of the compact operator G(H—\I)~ 1 G. 
Then for a = fi~ l , the corresponding eigenfunction it satisfies equation (|4.10[) . □ 

Now we can construct a counterexample to estimate (jl.ljl with arbitrary 5 < 
•v/a(A) for eigenvalues in gaps. 

Counterexample 4.4. Let a periodic function be such that operator ()4.1j) 

has infinite number of gaps. Then for an arbitrary e > 0, there exists an eigenvalue 
A of operator (|4. 5|) where Q S C^°(M) such that the corresponding eigenfunction 
equals 

${x) = e- slxl P± {x) (4.11) 

with ^ < E\J d(\) for sufficiently large ±x > 0. Here both p±(x) are either periodic 
or antiperiodic functions which are not identically zero. 

Indeed, using Proposition ^. 21 we can find A such that lnp(A) < £y/d(X). Then, 
by Lemma [4.31 we can construct Q(x) such that operator (|4.5p has eigenvalue A. 
The corresponding eigenfunction tp(x) is given by formula (|4.6[) for sufficiently large 
±x > 0. By virtue of (|4. 3[) this leads to representation (|4.11[) with 5 = lnp(A). 

It is possible that the results of this section can be deduced from very general 
results of [5]- 

Thanks are due to N. Firsova for a useful discussion of the periodic problem. 
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